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1 Abstract 

An exact solution is obtained for the electromagnetic field due to an electric current in the presence 
of a surface conductivity model of graphene. The graphene is represented by an infinitesimally- 
thin, local and isotropic two-sided conductivity surface. The field is obtained in terms of dyadic 
Green's functions represented as Sommcrfeld integrals. The solution of plane-wave reflection and 
transmission is presented, and surface wave propagation along graphene is studied via the poles 
of the Sommerfeld integrals. For isolated graphene characterized by complex surface conductivity 
it = a' +ja", a proper transverse-electric (TE) surface wave exists if and only if a" > (associated 
with intcrband conductivity) , and a proper transverse-magnetic (TM) surface wave exists for a" < 
(associated with intraband conductivity). By tuning the chemical potential at infrared frequencies, 
the sign of a" can be varied, allowing for some control over surface wave properties. 



2 Introduction 



Fundamental properties and potential applications of carbon-based structures are of interest in 
emerging nanoelectronic applications. Particularly promising is graphene, which is a planar atomic 
layer of carbon atoms bonded in a hexagonal structure. Graphene is the two-dimensional version 
of graphite, and a single- wall carbon nanotube can be thought of as graphene rolled into a tube pQ . 

Only recently has it become possible to fabricate ultrathin graphite, consisting of only a few 
graphene layers [2]-[3], and actual graphene [l]-[8]. In graphene, the energy- momentum relation- 
ship for electrons is linear over a wide range of energies, rather then quadratic, so that electrons 
in graphene behave as massless relativistic particles (Dirac fermions) with an energy-independent 
velocity. Graphene's band structure, together with it's extreme thinness, leads to a pronounced 
electric field effect [5], [9], which is the variation of a material's carrier concentration with electro- 
static gating. This is the governing principle behind traditional semiconductor device operation, 
and therefore this effect in graphene is particularly promising for the development of ultrathin car- 
bon nanoelectronic devices. Although the electric field effect also occurs in atomically thin metal 
films, these tend to be thermodynamically unstable, and do not form continuous layers with good 
transport properties. In contrast, graphene is stable, and, like its cylindrical carbon nanotube ver- 
sions, can exhibit ballistic transport over at least submicron distances [5]- Furthermore, it has been 
shown that graphene's conductance has a minimum, non-zero value associated with the conductance 
quantum, even when charge carrier concentrations vanish [6]. 

In this work, the interaction of an electromagnetic current source and graphene is considered. 
The electromagnetic fields are governed by Maxwell's equations, and the graphene is represented by 
a conductivity surface 10J that must arise from a microscopic quantum-dynamical model, or from 
measurement. The method assumes laterally infinite graphene residing at the interface between 
two dielectrics, in which case classical Maxwell's equations are solved exactly for an arbitrary elec- 
trical current. Related phenomena are discussed, such as plane-wave reflection and transmission 
through graphene [11] , and surface wave excitation and guidance, which is relevant to high-frequency 
electronic applications. It is found that the relative importance of the interband and intraband con- 
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tributions to the conductivity dictate surface wave behavior [12) , and that surface wave propagation 
can be controlled by varying the chemical potential. 

Although at this time only graphenc samples with lateral dimensions on the order of tens of 
microns have been fabricated, the infinite sheet model provides a first step in analyzing electro- 
magnetic properties of graphenc. It is also relevant to sufficiently large finite-sized sheets, assuming 
that electronic edge effects and electromagnetic edge diffraction can be ignored. In the following 
all units are in the SI system, and the time variation (suppressed) is e 3Ult , where j is the imaginary 
unit. 



3 Formulation of the Model 

3.1 Electronic Model of Graphene 

Figure 1 depicts laterally-infinite graphene lying in the x — z plane at the interface between two 
different mediums characterized by /ii,£i for y > and 1x2^2 for y < 0, where all material 
parameters may be complex- valued. 

The graphenc is modeled as an infinitesimally-thin, local two-sided surface characterized by 
a surface conductivity a (oj, /i c ,r , T), where u> is radian frequency, fi c is chemical potential, T 
is a phenomenological scattering rate that is assumed to be independent of energy s, and T is 
temperature. The conductivity of graphene has been considered in several recent works [llj— |18j. 
and here we use the expression resulting from the Kubo formula [19) , 



a [u, He, T, T) = — 



1 r 4m- a -^^ a) 



(oj - j2T) 2 Jo V de de 



fd (-£) - fd (e) ± 
{LO-J2T) 2 -A{e/Hf 



where e is the charge of an electron, H = h/2n is the reduced Planck's constant, fd (e) — {e^ E ^ tJ -^/ kBT + l) 
is the Fermi-Dirac distribution, and ks is Boltzmann's constant. We assume that no external mag- 
netic field is present, and so the local conductivity is isotropic (i.e., there is no Hall conductivity). 
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The first term in (JTJ) is due to intraband contributions, and the second term to interband contribu- 
tions. 

For an isolated graphene sheet the chemical potential /i c is determined by the carrier density 

n„, 
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3C 



n s = -£n2 e (/d (e) - /d (e + 2/ic)) de, (2) 

TTIl-Vp Jq 

where vf — 9.5 x 10 5 m/s is the Fermi velocity. The carrier density can be controlled by application 
of a gate voltage and/or chemical doping. For the undoped, ungated case at T — K, n s = fi c = 0. 
The intraband term in ([1]) can be evaluated as 

"»■« <«• F - r » - -i^m (3) 

For the case = 0, (J3j) was first derived in [20] for graphite (with the addition of a factor to 
account for the interlayer separation between graphene planes), and corresponds to the intraband 
conductivity of a single-wall carbon nanotube in the limit of infinite radius [10] . With a = a' +ja", 
it can be seen that <j' intra > and cr'/ ntra < 0. As will be discussed later, the imaginary part of 
conductivity plays an important role in the propagation of surface waves guided by the graphene 
sheet [12] . 

The interband conductivity can be approximated for kgT <§C |/x c | , fko as [2Tj . 

a mter ( W , ^ F, 0) * — In [ 2Wl + {uj _ j2r)n ) , (4) 

such that for r = and 2 |/x c | > fej, (7i„ ter = jcr'/ nter with cr-'„ ter > 0. For r = and 2 \fi c \ < fkj, 
Winter is complex-valued, with [52] 



. ,„ 1U = 6.085 x 10~ 5 (S), (5) 



and cr" nter > for fi c ^ 0. 
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3.2 Dyadic Green's Function for a Surface Model of Graphene 

For any planarly layered, piecewise-constant medium, the electric and magnetic fields in region n 
due to an electric current can be obtained as [23]-[24] 

BW (r) = (fc 2 + VV-) ^ (n) (r) , (6) 
(r) = jue n V x (r) , (7) 

where k n — u>^p n e n and ir^ (r) are the wavenumber and electric Hertzian potential in region n, 
respectively. Assuming that the current source is in region 1, then 



^ 1 )(r)=^(r)+ 7 r 1 (r)= / {g? (r,r') + g* (r,r')| • dtf, 



(8) 



^( 2 > (r) = «> (r) = / g* (r, r') • dfi', (9) 

where the underscore indicates a dyadic quantity, and where f2 is the support of the current. With 
y parallel to the interface normal, the principle Greens dyadic can be written as |23j 



p-jkiR i roo u-(2) /■ s 

g P(r ) r')=I- =1—1 e -^\v-y'\ th SMI kdk (i ) 



where 



pi = K-K, P = ^{x-x'Y + {z-z'y, (ii) 



R=\v-v'\ = ^(y-y>y+p*, 

and where fc p is a radial wavenumber and I is the unit dyadic. 

The scattered Green's dyadics can be obtained by enforcing the usual electromagnetic boundary 
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conditions 



yx (H 1 -H 2 )=J«, 

yx(E 1 -E 2 )=-j;, (12) 

where J* (A/m) and (V/m) are electric and magnetic surface currents on the boundary. For a 

local model of graphene in the absence of a magnetic field and associated Hall effect conductivity, 
a is a scalar. Therefore [TO] , 

J tx ( x > y = 0,z) = aE x (x, y = 0,z), (13) 

J s e ,z 0; y = 0,z) = <tE z (x, y = 0,z), (14) 

3 s m (x,y = 0,z) = 0, (15) 

and ([T2]) becomes 

(y = 0+) = E 2 , a (y = 0") , a = x,z, (16) 

H 2 , x (y = 0") - (y = 0+) = aE z (y = 0) , (17) 

#2, z (» = 0") - ff M (y = 0+) = -oE x {y = 0) . (18) 

Using ([6]), the boundary conditions on the Hertzian potential at (x,y = 0, z) are 

n ha = 7V 2 M 2 ^ 2iQ , (19) 

(7 

eiTTl.y - £2^2,y = — V • 7Tl (20) 

£ 1^7^ = —kt^i.a (21) 

dy dy 

a = z, where iV 2 = £2/^1 and M 2 = fa/vi- I n the absence of magnetic contrast (e.g., if M = 1) 
and surface conductivity, boundary conditions l|19 p - (j22j) are identical to the Hertzian potential 

G 



boundary conditions presented 25J. 

Enforcing (fT9)) - (|2"2")) and following the method described in [35], the scattered Green's dyadic is 
found to be 



g* (r, r') = yy g s n (r, r') + ( yx — + yz— ) g s c (r, r') + (xx + zz) g s t (r, r') , (23) 



where the Sommerfeld integrals are 



2tt 



h ^ (h n\ P -pi{y+y') 
4pi 



(24) 



(3 = t,n, c, with 



Rt 

Rn 

Rr 



M 2 p 1 ~p 2 ~ iouj[i 2 _ N H (k p ,u>) 
M 2 pi +p 2 + j<Jujfi 2 Z H (k p , uj) ' 

N E {k p ,w) 
Z E (k 



N 2 pi - p 2 



2pi 



N l M 2 - 1 + 



ap 2 M 2 



Z H Z E 



(25) 
(26) 

(27) 



which reduce to the previously known results as a — > 0. 

The Green's dyadic for region 2, g* (r, r'), has the same form as for region 1, although in (I24| 
the replacement 

R f3 e- pi ( y+y ') -> T/3 e P2V e- piv ' (28) 



must be made, where 



T t = 
T n = 



(l + Rt) 2p x 



N 2 M 2 N 2 Z H ' 
Pi (1 - Rn) _ 2P1 



I>2 



2pi 



(N 2 M 2 - 1) + 4^ 

\ ) JUS! 



N 2 Z H Z E 



(29) 
(30) 

(31) 
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As in the case of a simple dielectric interface, the denominators Z H ' E (k p , u>) = implicate 
pole singularities in the spectral plane associated with surface wave phenomena. Furthermore, 



both waveparameters p n = yk p — fc^, n = 1,2, lead to branch points at k p — ±fc n , and thus 
the fcp-plane is a four-sheeted Riemann surface. The standard hyperbolic branch cuts [23) that 
separate the one proper sheet (where Re(p„) > 0, such that the radiation condition as \y\ — ► oo 
is satisfied) and the three improper sheets (where Re(p„) < 0) are the same as in the absence of 
surface conductivity a. 

In addition to representing the exact field from a given current, several interesting electromag- 
netic aspects of graphene can be obtained from the above relations. 

3.3 Plane- Wave Reflection and Transmission Coefficients 

Normal incidence plane-wave reflection and transmission coefficients can be obtained from the 
previous formulation by setting k p = in (|2"5)) and (|2l)|) . To see this, consider the current 

J«(r)=S^°5(r-ro), (32) 

UfJ.1 

where a = x or y, ro = yyo, and where j/o 3> 0. This current leads to a unit-amplitude, S-polarized, 
transverse electromagnetic plane wave, normally-incident on the interface. The far scattered field 
in region 1, which is the reflected field, can be obtained by evaluating the spectral integral (|24p 
using the method of steepest descents, which leads to the reflected field E r = aTe~i lV , where the 
reflection coefficient is T — R t (k p = 0). In a similar manner, the far scattered field in region 2, the 
transmitted field, is obtained as E* = aTe^ 211 , where the transmission coefficient is T — (1 + T). 
Therefore, 



T= m-m-°rnr» (33) 

T=(l + R t )= r (34) 

[V2 + m + vvim) 
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where r\ n = y/ \i n je n is the wave impedance in region n. The plane-wave reflection and transmission 
coefficients obviously reduce to the correct results for a = 0, and in the limit a — > oo, T — > — 1 and 
T — > as expected. 

In the special case ei = £2 = £0 an d /zi = p,2 = MOj 



"in 
2 



where 770 = \/ Mo / to — 377 ohms. The reflection coefficient agrees with the result presented in [11] 
for normal incidence. 

3.4 Surface Waves Guided by Graphene 

Pole singularities in the Sommerfeld integrals represent discrete surface waves guided by the medium 
[23]-[2l]. From (|25 |) -([27 |) and (f29|) -(l3T |) . the dispersion equation for surface waves that are transverse- 
electric (TE) to the propagation direction p (also known as -ff -waves) is 

Z H (k p , w) = M 2 Pl + P2 + jauj m = 0, (36) 

whereas for transverse- magnetic (TM) waves (E- waves), 

Z E (k p , w) = N 2 Pl + P2 + ^™ = 0. (37) 

In the limit that £1 — £2 — £0, Z H ' E agree with TE and TM dispersion equations in [T2] 

The surface wave field can be obtained from the residue contribution of the Sommerfeld integrals. 
For example, the electric field in region 1 associated with the surface wave excited by a Hertzian 
dipole current J (r) = yA 5 (x) 5 (y) S (z) is 

(i) f ^ Aok P R 'n \ , e±\ ttW ru ^ Jkl+pl)H^ (k pPo ) \ 



E (po) = —^e'^y x-+i- H^' (k pP0 ) y 
4w£i V po po 1 



H (kppo) - y V , 

/ kjkl-kj J 
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where R' n = N E / (dZ E /dk p ) , H { 2) ' (a) = dH {2) (a) /da, and p = Vx 2 + z 2 . The term e - piy leads 
to exponential decay away from the graphene surface on the proper sheet (Re(p n ) > 0,n = 1,2). 
The surface wave mode may or may not lie on the proper Riemann sheet, depending on the value 
of surface conductivity, as described below. In general, only modes on the proper sheet directly 
result in physical wave phenomena, although leaky modes on the improper sheet can be used to 
approximate parts of the spectrum in restricted spatial regions, and to explain certain radiation 
phenomena [26] , 

3.4.1 Transverse-Electric Surface Waves 

Noting that p 2 — p\ = k 2 , {n\e\ — ^2^2): where [i r n and e r n are the relative material parameters (i.e., 
/!„ = iJL r n Ho and e n — e r n eo) and /cq = w 2 /j,o£o is the free-space wavenumber, then if M = 1 (fi[ — 
/i 2 = Mr) t ne TE dispersion equation (|36|) can be solved for the radial surface wave propagation 
constant, yielding 

k ? = fc ° V M " £l - { 2^ J • (38) 

If, furthermore, N — 1 [e\ = e T 2 = e r ), then (|38|) reduces to 

k p = fcoyMr£r - (^^p-^j . (39) 

For the case 1, then (|36|) leads to 

k P = fco^i - J^-^I (^ 2 ^oM2 T \/ (<w£) 2 - (M 4 - 1) (eM - eM)) • (40) 
Considering the special case of graphene in free-space, setting e\ = e r 2 = [i\ = [i\ I = 1, 



*P = *»V 1 -( f r)- w 

If ex is real-valued (cr = er') and (cr'770/2) 2 < 1, then a fast propagating mode exists, and if 
(cr'?7o/2) > 1 the wave is either purely attenuating or growing in the radial direction. How- 
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ever, in both cases p n = Po = y {k p ) 2 — kfi = —ja'ujfiQ/2 from and so Re (po) > is violated. 
Therefore, for isolated graphcne with a real- valued (i.e., when a = cr m i n , at low temperatures and 
small jjL c ), all TE modes are on the improper Riemann sheet. The fast leaky mode may play a role 
in radiation from the structure. 

If the conductivity is pure imaginary, a = jcr", then k p > ko and a slow wave exists. In this case 
Po = (r"ujfj,o/2 and if a" > then Re (pa) > and the wave is a slow surface wave on the proper 
sheet. This will occur when the interband conductivity dominates over the intraband contribution, 
as described in |12j . However, if a" < 0, which occurs when the intraband contribution dominates, 
the mode is exponentially growing in the vertical direction and is a leaky wave on the improper 
sheet. 

More generally, for complex conductivity, 

Po = ^ = ( CT ~ 3a >~' ( 42 > 

and therefore if a" < the mode is on the improper sheet, whereas if a" > a surface wave on the 
proper sheet is obtained. 

3.4.2 Transverse-Magnetic Surface Waves 

TM waves are governed by the dispersion equation (|37[) . For general material parameters this 
relation is more complicated then for the TE case, and so here we concentrate on an isolated 
graphcne surface (e\ = e r 2 = [i\ — = 1). Then, pq = —j2ujeo/<J and 



*, = fayi-(^)\ (43) 

If a is real-valued then Re (po) > is violated and TM modes are on the improper Riemann sheet. 

If conductivity is pure imaginary, then k p > ko and a slow wave exists. Since po = —2ujeo/a", 
if a" > then the wave is on the improper sheet, and if a" < the mode is a slow surface wave on 
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the proper sheet. For complex conductivity, 



pa = = 2 (o- + jo- ) , (44) 



and therefore if a" < (intraband conductivity dominates) the mode is a surface wave on the 
proper sheet, whereas if a" > (interband conductivity dominates) the mode is on the improper 
sheet. 

In summary, for isolated graphene a proper TE surface wave exists if a" > 0, resulting in the 
radial wavenumber (|4ip . and a proper TM surface wave with wavenumber ()43|) is obtained for 
a" < 0. The same conclusions hold for graphene in a homogeneous dielectric. For graphene with 
a" = 0, no surface-wave propagation is possible. Note that this only refers to wave-propagation 
effects; dc or low-frequency transport between electrodes can occur, leading to electronic device 
possibilities. In fact, in devices not based on wave phenomena, the absence of surface waves is 
usually beneficial, as spurious radiation and coupling effects, and the associated degradation of 
device performance, are often associated with surface wave excitation. 

The degree of confinement of the surface wave to the graphene layer can be gauged by defining 
an attenuation length £, at which point the wave decays to 1/e of its value on the surface. For 
graphene embedded in a homogeneous medium characterized by e and /i, C, — Re (p), leading to 
( TE = 2ja"uj[i (a" > 0) and C™ = - |cr| 2 /2ujea" (a" < 0). When normalized to wavelength, 

/-TE -i 

a" > (45) 



A Trrja" ' 

™ Ti\a\ 2 

n] 1 , a"<0. (46) 



A 471-cr" ' 

Obviously, strong confinement arises from large imaginary conductivity, as would be expected 
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4 Results 



In this section, some results are shown for surface wave characteristics of graphene in the microwave 
and infrared regimes. In all cases T = 0.11 meV, T = 300 K, and an isolated graphene surface (i.e., 
when the surrounding medium is vacuum) is considered. The value of the scattering rate is chosen 
to be approximately the same as for electron-acoustic phonon interactions in single-wall carbon 
nanotubes [27] . 

We first consider the case of zero chemical potential at microwave and far-infrared wave frequen- 
cies. Fig. 2 shows the complex conductivity, TM surface- wave wavenumber (|43[) . and attenuation 
length (|46[) . In this case the intraband conductivity is dominant over the interband contribution, 
and so a" < 0, such that only a TM surface wave can exist. The dispersion of the complex conduc- 
tivity follows simply from the Drude form (|3|). At low frequencies the TM surface wave is poorly 
confined to the graphene surface (C™/A ^ 1); an d therefore it is lightly damped and relatively 
fast (i.e., fcJ M /fco ~ 1). As frequency increases into the far-infrared, the surface wave becomes more 
tightly confined to the graphene layer, but becomes slow as energy is concentrated on the graphene 
surface. 

The conductivity can be varied by adjusting the chemical potential, which is governed by the 
carrier density via The carrier density can be changed by either chemical doping or by the 
application of a bias voltage. Fig. 3 shows the interband and intraband conductivity as a function 
of chemical potential for ui — 6.58 /LteV (ui = 10 GHz). As expected, the conductivity increases 
with increasing chemical potential, associated with a higher carrier density n s . Because intraband 
conductivity is dominant, a" remains negative as chemical potential is varied, and therefore only 
a TM surface wave may propagate. The insert shows the real part of the intraband conductivity 
(the dominant term) on a linear scale, showing the linear dependence of a on \i c . 

Since the TM surface-wave is not well confined to the graphene surface at lower microwave 
frequencies, it undergoes little dispersion with respect to chemical potential. This is illustrated in 
Fig. 4, where the TM wavenumber and attenuation length are shown as a function of chemical 
potential at u> = 6.58 ^eV. Because of the simple form for the TM surface- wave wavenumber (|4"3"1) . 
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the wavenumber and attenuation length merely follow the conductivity profile. 

However, at infrared frequencies moderate changes in the chemical potential can significantly 
alter graphene's conductivity, and, significantly, change the sign of its imaginary part. Fig. 5 shows 
the various components of the graphene conductivity at u> = 0.263 eV (w = 400 THz), and in Fig. 
6 the total conductivity (interband plus intraband) is shown. From ([¥]), at T = T = an abrupt 
change in Omter occurs when 2 \/j, c \ = huj, which in this case is \/x c \ = 0.132 eV, denoted by the 
vertical dashed lines in the figures. Since the associated Fermi temperature is several thousand K, 
the T — behavior qualitatively remains the same at 300 K, although the discontinuity is softened 
due to the higher temperature. It can be seen that for |/z c | less than approximately 0.13 eV, <j" nter 
dominates over cr" ntra and a" > 0, so that only a proper TE surface wave mode exists. Outside 
of this range, only a TM surface wave propagates. This is shown in Fig. 7, where it can be seen 
that the TM mode is moderately dispersive with chemical potential, especially in the region of the 
sign change in a" near |/i c | ~ 0.132 eV, since the TM surface wave is fairly well confined to the 
graphene surface ((™ /\ < 10~ 2 ). In the region where a" is positive, the TE mode exists but 
is poorly confined to the graphene surface, and so it is essentially nondispersive and very lightly 
damped. The oscillatory behavior of the attenuation length C, TE follows simply from the form of 
a" via 05]). 

In Fig. 8 the conductivity is shown as a function of frequency in the infrared regime for a fixed 
value of chemical potential, \x c — 0.1 eV, at 300 K. The dispersion behavior of the conductivity 
follows simply from ([3]) and ((4]). The point 2 |/x c | = fko occurs at u) = 0.2 eV (u) ~ 301.6 THz), 
whereupon the interband contribution dominates and a" becomes positive (the intraband contri- 
bution varies as uj^ 1 for uj T, and so becomes small at sufficiently high infrared frequencies). For 
comparison, the T = result is also shown. 

Fig. 9 shows the TM and TE wavenumbers and attenuation lengths for the conductivity profile 
given in Fig. 8. For w < 0.2 eV only a TM surface wave exists since a" < 0, and for ui > 0.2 
eV only a TE surface wave exists due to a" > 0. From Fig. 8 it is clear that above u> = 0.1 eV, 
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\a\ 770/2 <C 1, and so 




The TM mode is tightly confined to the graphene surface, and from (|4"T)) it can be seen that in 
vicinity of the transition point at uo = 0.2 eV, where a" ~ 0, k 1 ' will be large and the TM mode will be 
highly damped. Just to the right of the transition, from (1481) the TE wavenumber is predominately 
real, and so the TE mode is very lightly damped. From (|4"S")) it is clear that ( TE /A 3> 1, and so the 
TE mode is not well confined to the surface in this frequency range. 

5 Conclusions 

An exact solution has been obtained for the electromagnetic field due to an electric current near a 
surface conductivity model of graphene. Dyadic Green's functions have been presented in terms of 
Sommerfeld integrals, plane-wave reflection and transmission coefficients have been provided, and 
surface-wave propagation on graphene has been discussed in the microwave and infrared regimes. 
The relative importance of interband and intraband contributions for surface wave propagation has 
been emphasized. 
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Figure Captions: 

Fig. 1. (a) Depiction of graphcne (top view), where the small circles denote carbon atoms, and 
(b) graphene characterized by conductance a at the interface between two dielectrics (side view). 

Fig. 2. Complex conductivity, TM surface-wave wavenumber, and attenuation length for fi c — 
at 300 K in the microwave through far-infrared frequency range (ui = 15.2 GHz to 15.2 THz). 

Fig. 3. Intcrband and intraband conductivity as a function of chemical potential at T = 300 
K, u) = 6.58 [icV (oj = 10 GHz). Note the logarithmic scale; the insert shows the real part of the 
intraband conductivity on a linear scale, showing the linear dependence of a on /i c . 

Fig. 4. Attenuation length and surface-wave wavenumber for the TM mode as a function of 
chemical potential at T = 300 K, uj — 10 GHz (6.58 /xeV). The corresponding conductivity profile 
is shown in Fig. 3. 

Fig. 5. Interband and intraband conductivity as a function of chemical potential at T — 300 K, 
oj = 0.263 cV (w = 400 THz). The dashed vertical lines represent the point 2 \fi c \ — ftu. 

Fig. 6. Total (interband plus intraband) conductivity as a function of chemical potential at 
T = 300 K,uj = 0.263 eV (w = 400 THz). The dashed vertical lines represent the point 2 \fj, c \ = Hu. 

Fig. 7. Attenuation length and surface-wave wavenumbers as a function of chemical potential 
at T = 300 K, ui = 0.263 cV (ui = 400 THz). The corresponding conductivity profile is shown in 
Figures 5 and 6. TE and TM modes are shown, although only portions where wavenumbers lie on 
the proper Riemann sheet are provided. 

Fig. 8. Total conductivity (interband plus intraband) as a function of frequency at T — 300 K 
and fi c = 0.1 eV at infrared frequencies (uu = 0.1 cV~ 151.2 THz). The T = result is also shown 
for comparison. 

Fig. 9. Attenuation length and surface-wave wavenumbers as a function of frequency at [i c = 0.1 
eV, T = 300 K. The corresponding conductivity profile is shown in Fig. 8. 
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